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Abstract 


A simple model of solar flare radiation event 
particle transport is developed to permit the cal- 
culation of fluxes and related quantities as a 
function of distance from the sun (r). This model 
assumes the particles spiral around the solar mag- 
netic field lines with a constant pitch angle. 

The particle angular distributions and onset plus 
arrival times as functions of energy at 1 AU agree 
with observations if the pitch angle distribution 
peaks near 90^. g consequence the time depend- 
ence factor is essentially proportional to R"1*T^ 
(R in AU), and the event flux is proportional to^ 
R“^. 


O' ) 

n I approximately constants 
m \ 


If the characteristics of the largest solar flare 
particle events are used to evaluate a , n, and 
m, the values obtained are approximately 0.022, 
O.k, and 0.75, respectively. An attempt to 
incorporate the variations of these quantities 
due to radiation event size led to the relation- 
ships'^ ^ ' 
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I. Introduction 


Solar flare particle events constitute one 
of the important sources of natural nuclear radia- 
tion at 1 AU from the sun. Approximately 100 such 
events have been observed at the earth (see Table 
ranging up to - lO^^ protons /cm^.event 


above 10 Mev. Various particle flux models have 
been developed describing the time and energy 
behavior observed at 1 AU. These models, coupled 
with prediction mechanisms (usually correlated to 
sunspot numbers) have been used to estimate the 
expected solar flare particle environment at 1 AU 
for future space missions. These techniques have 
also been used for missions to Mars (I.5 AU) and 
Venus (0.7 AU). 


For missions which involve sending spacecraft 
far from 1 AU (e.g., the Mercury- Venus and the 
Outer Planet missions), the solar flare particle 
model developed from observations at 1 AU is not 
adequate. Specifically, the dependence of the 
solar flare particle event model on distance from 
the sun must be incorporated. This paper presents 
a simple model for incorporating this spatial 
dependence which agrees with the observations at 
1 AU. 


II. Model at 1 AU 

The flux of solar flare p>articles observed at 
1 AU is observed to rise qua si -linearly over a 
period of some hours to a peak value, then decay 
approximately exponentially over an appreciably 
longer period. The onset, rise, and characteris- 
tic decay times are functions of particle energy. 

A previously developed mathematical e:«)ression 
which has these characteristics 

t) = Ate-^Ct 

cm^-hr 

where ig the particle energy(Mev) 

t is the time from the arrival of the first 
particles (hrs) 

A is a normalization factor 


where all quantities are as previously defined. 
Incorporation of these relationships in Equation 
(1) yields formulae which account for the shorter 
rise and decay times and steeper energy spectra 
for the smaller events. However, the smaller 
event characteristics vary so much that generali- 
zations are difficult and attempts to account for 
the size-dependent characteristics are generally 
more trouble than they are worth. 

Based on Equation (l), the following relations 
have been derived: 
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These relationships have proven useful when 
applied to large solar flare radiation events at 
1 AU. In particular, they have been applied to 
events for which the data was incomplete— some 
of the numbers in Table I were obtained this way. 

These formulae have no dependence on the dis- 
tance frcan the sun, and are therefore strictly 
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applicable to 1 AU only- In the following 
sections, modifications of these formulae to 
account for distance from the sun are discussed. 

III. Particle Transport 

In order to account for the dependence of 
the solar flare particle event model on distance 
from the sun, it is necessary to have a particle 
transport model. Since particle acceleration and 
particle transport are related during the accel- 
eration portion of the particle event, it is 
desirable to go back to the solar flare itself. 

A solar flare is an event observed to take 
place above the photosphere of the sun, generally 
lasting < 20 minutes. The flare, which is most 
readily observed at the Lynan - O' wavelength of 
ionized hydrogen (-1215 A) is believed to derive 
its energy from the collapse of the magnetic 
fields associated with sunspots. Various theories 
have been advanced to account for the acceleration 
of charged particles during or as a result of a 
solar flare. Among these theories are betatron 
acceleration and shock wave acceleration. In 
betatron acceleration, a charged particle orbiting 
above a collapsing sunspot gains energy as in a 
laboratory particle accelerator of the same name. 
In shock wave acceleration the absorption of the 
Lyman - a photons in the inner corona produces a 
supersonic blast wave which accelerates the 
charged particles. These proton acceleratipu 
mechanisms have been investigated by Gold,'^"^>^ 
Wentzel,(5) Parker,^ Weddell,^ and others. 

The evidence is not conclusive for any single 
mechanism, and both betatron and shock wave 
accelerations probably contribute. 

Once the charged particles have been accele- 
rated, their propagation away from the sun is 
controlled by the interplanetary magnetic fields. 
These interplanetary fields are complex, having 
time and spatial dependences which are not well 
known. However, by making some simplifying 
assumptions, it is possible to obtain a spatially 
dependent solar flare particle model. 

If it is assumed that the quiet sun inter- 
planetary magnetic field is essentially undis- 
turbed during the propagation of the solar flare 
protons, it is possible to write: 



r 


J = m Vj^ r 
= V sin P 


is the aximuthal component of B (webers/m^) 

R is the distance from the sun (AU) 
q is the charge of the particles (coulombs) 

V is the velocity of the charged particles 

(meters/sec) , 

V is the perpendicular component of v (meters/ 

^ sec ) 

v,| is the parallel component of v (meters/sec) 

p is the pitch angle (angle v makes with the 
guiding center line of B). Since ^ 

10"^ webers/m^ ( 3.5 gammas) at 1 AU, B is 5 
gammas, and the gryroradius is 2. 75 x 10 «^/e 
meters, where E is in Mev. The gyroradius as a 
function of distance from the siin is approximately 

r = 3.9 X lo"^ r E ^ meters (6) 

y 1 + r2 

This function is plotted in Figure 2. 

Assuming the particles spiral around the mag- 
netic lines of force, it is expected that their 
gyroradius will increase as their distance from 
the sun increases. Based upon the above equations 
it is expected that their pitch angle p will 
decrease with distance from the sun, according 
to the well-known relationship. 


J q sin^P ^ . 

= = constant 

2 m E B 

However, the pitch angle ( P) cannot exceed 90^. 

If the charged particle has this limiting angle 
close to the sun, its pitch angle might be 
expected to become quite small by the time it 
reaches 1 AU. If a particle has a pitch angle 
of 90 ° starting at 10 solar radii, the calculated 
pitch angle will be — 3° by the time it reaches 
1 AU (-^215 solar radii) according to this model. 
Since solar flare event particles are observed 
at 1 AU to be essentially isotropic after the fl\^ 
peak, this simplified model is clearly unrealistic. 

It is possible to obtain reasonable agreement 
with observations at 1 AU by making the assumption 
that the pitch angle is independent of distance 
from the sun. This accounts for the fact that 
particles of a given energy with small pitch 
angles will arrive first (as is observed), and 
the resultant particle angular distributions 
spread as functions of time. To the extent that 
this model appears to agree with the observations 
at 1 AU, it may be useful for calculating the 
spatial dependence of solar flare particle radia- 
tion. 

The transit time (t^^) as a function of parti- 
cle energy and pitch angle (both assumed to be 
constants of the motion) is 


Vjj = V cos p 
where 

B is the solar magnetic field (webers/m^) 

B is the radial component of B (webers/m ) 
^ (see Figure l) 



where s is the length of the line of magnetic flux 
which serves as the guiding center for the i)arti- 
cle gyrorotation. This length(s) maybe calcu- 
lated from the relationship: 
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But since 


=y^ds =J^{Tr)2 + (r de )' 


tan“^ r de _ +an“^ Ti j. -1 

- xan Bq tan r 


r d9 = R dr 


Therefore 


s = I dr. /I + R" 


= - j^iyi . R + log (R yi ^ r2^ j 

When the distance from the sun (r) is small the 
two terms are approximately equal in value, but 
as the distance from the sun Increases, the first 
term predominates. This function is graphically 
displayed in Figure 3. A reasonable approxima- 
tion is: 



If it is assumed that the pitch angle (p) of the 
particle is a constant of the motion, the transit 
time as a function of pitch angle and particle 
energy may be calculated. Based on Equation (7), 
the resxilt is (for protons): 

^tr 5-5 X 10^ [r a + R^'+ log (r + >/i + R^ )j 

VF cosP 

seconds 

where R is in AU and E is in Mev. With only a 
few percent error in the region of interest (O.l - 
10 AU) this may he approximated by 


1.1 x.lO^ 


y/W* cos p > 


seconds 


Values of this function are listed in Tables II- 
IV, and are plotted in Figures (4) to (6). These 
figures show the expected dependence on proton 
energy (e), pitch angle ( P ), arid distance from 
the sun (R). In principle, this expression may 
be solved to obtain any of the quantities as 
functions of the others. The quantity of major 
interest is the proton flux as a function of these 
parameters. This is determined by the initial 
proton distribution in p , assuming no protons 
are lost after being accelerated. 

If no particles are lost after being accele- 
rated, the total number per cm^ per event is 
expected to be proportional to R“2 ^ distance 

from the sun in AU). Assuming an r 2 time-inte- 
grated spatial dependence effectively makes the 
solar flare particle event probability spatially 
independent. Any other assumption introduces the 
complication of a spatially dependent probability 


function. Since the time- integrated particle flux 
for a mission is essentially the product of the 
flux per event and the probability of that event, 
little, if anything, is gained by considering a 
spatially dependent particle event probability. 

If a significant fraction of solar flare 
particles is lost after being accelerated, the 
flux per event will exhibit an R"^ spatial 
dependence where n*> 2 over the region where 
particle loss is important. The major probable 
loss mechanism is collision with objects in space 
(magnetospheres, planets, moons, etc.) Consider- 
ing the relative emptiness of space this is not 
expected to be a major loss mechanism. Even the 
asteroid belt is only^lO-6 

that the solar flare particle traversing it will 
be small, (it may not be completely negligible, 
however, since the proton helical path length 
through the asteroid belts will be - two orders 
of magnitude larger than that of a solar photon) . 
The probability of a solar flare particle encoun- 
tering a planet, even assuming the particles are 
essentially confined to the ecliptic plane is<l^. 

It is interesting to calculate the pitch angle 
which corresponds to the peak flux rate. At 1 AU 
the onset + rise time (time from the flare on the 
sun to the peak flux rate) is approximately: 

^rise ^ hrs 

E OA 
o 

where Eq is the proton energy in Mev. (This 
relationship was obtained by averaging the values 
observed for the largest radiation events— smaller 
events often have shorter onset + rise times). 

Since the proton transit time (given by Equation 
9) is proportional to E-0-5^ there is apparently 
a constant delay time to be added. A little 
arithmetic shows that this delay time is 86o/cos 6 
seconds. Thus: 

t. . = + • 860 

"rise ^'tr seconds 

cos p 

This accounts for the acceleration time during 
which the protons do not migrate away from the sun 
significantly as well as the time between the 
arrival of the first protons of a given energy 
and the peak flux of these protons. It will be 
noted that a fairly steep proton energy spectriun 
has been assumed (based on observations at 1 AU), 
since the onset + rise time for particles above 
energy E^ is taken to be essentially that of 
particles of energy E^, 

For 1, 10, and 100 Mev protons at 1 AU, the 
observed onset + rise times are ~ 1.6 x, 

10 seconds, ^6.5 x 10^ seconds, and~ 2.6 x 10^ 
seconds, respectively. The corresponding transit 
times are ~ 1 5 x lo5 seconds, ~5.2 x 10»^ seconds, 
and <-1.5 x 10^ seconds. These are the transit 
times for protons with pitch angles of - 86°. it 
will be noted that any change of pitch angle as 
the protons migrate away from the sun will decrease 
the transit time below that observed. Thus, unless 
there is some mechanism which acts to increase the 
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pitch angles of protons as they move away from 
the sun, the assumption of constant pitch angle 
appears to he superior to any obvious simple 
alternative. Also noteworthy, is that beyond 
2 AU solar flares on the back side of the sun can 
be important sources of particles. 

Another interesting conclusion, valid to the 
extent that the assumption of a constant pitch 
angle is valid is that the bulk of the protons 
have pitch angles close to 90°. Based upon 
Equation (2), only -26fo of the protons above a 
given energy have arrived at 1 AU by the time the 
peak flux rate of these protons has been reached. 
Thus , -- of the protons observed at 1 AU 
appear to have pitch angles between 86° and 90°. 
Assuming a random distribution of proton velocity 
directions prior to acceleration, a sin P pitch 
angle distribution would be expected from solid 
angle considerations. Apparently, the accelera- 
tion mechanism favors particles with pitch angles 
close to 90°, This seems to favor the betatron 
acceleration mechanism over the shock wave accel- 
eration mechanism. However, since the constant 
pitch angle assumption is merely a simplified 
model, any conclusions based upon it must be con- 
sidered uncertain. 
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The maximum pitch angle determines the limits 
of the particle angular distril^ution as a function 
of time and energy. The particle angular dis- 
tribution within these time and energy dependent 
limits according to this model is determined by 
the initial particle pitch angle distribution. 
However, not enough angular distribution infor- 
mation is available to derive statistically mean- 
ingful pitch angle distributions as functions of 
time and energy. When such data becomes available, 
its incorporation into this model is straight- 
forward. 


IV. Spatially Dependent Model 

It is now xxDSSible to use the transport 
model to calculate the parameters of the solar 
flare particle radiation as functions of distance 
from the sun. Since the pitch angle is assumed 
to remain constant, so does Vj| (the particle 
velocity parallel to the line of magnetic force 
which serves as the gyrorotation guiding center). 
The time for a particle to reach a given distance 
from the sun is simply proportional to the spiral 
path length (given by Equation ?)• Thus, the 
time dependence of the solar flare particle flux 
is expected to be 


t»-- 


f (R) 


where f (r) is given by Equation (?) and t^ is - 
1.2 X 10^ sec. It is, therefore, possible to 
rewrite Equations (l) to (5) incorporating R (the 
distance from the sun, as a parameter). The 
results are:^°i m -1 
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A (t-to) e 


R f"^(R) Eg“ 
-0.022Eq°-'^ (R)t 


R^ f^(R) E 
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particles 

cm^-hr 


-0.022E°-^f‘^ (R)t 

4>(E, t) = At e 

r2 

fo.T5 + 0.0088tf^(R)E°‘^] particles (1T) 

^ cm^-hr-Mev 

where all symbols have been -previously identified. 
As mentioned previously, f (R) can be approxi- 
mated by 

f (R)~R^l + ^ 

l.T 

For R >1, this may be further approximated by R 
with < 15^ error at 1 AU and a smaller error for 
KRO* 

The relative spatial dependences of the time- 
integrated particle flux rate ( « R f (R) 
are shown in Figure 7* These curves show that the 
relative solar flare particle environment becomes 
severe as the sun is approached. For the largest 
such particle event observed at 1 AU ( 11-12-60), 
the peak and time -integrated proton fluxes above 
10 Mev at Mercury (~0.^ AU) would be— 7*2 x 10 
p/cm^-sec and -6.3 x 10^^ p/cm^, respectively — 
enough to affect sensitive spacecraft components. 
For this same event the corresponding numbers at 
Jupiter C'^5*2 AU) would be 120 p/cm^-sec and 
3.7 X 10° p/cm^, which should cause no trouble. 

It will be noted that particle loss and- 
change of particle energy after acceleration have 
been neglected. If either of these changes takes 
place (and there is some evidence that they do), 
the solar flare particle environment as the .sun 
is approached will be even more severe than here 
calculated. Conversely, at distances greater than 
1 AU, the environment may very well be less severe 
than expected. 
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V. Conclusions 


Parameter in solar flare particle event 
model 


n 


The solar flare particle event model developed 
to fit the observed data at 1 AU has been modified 
to account for the distance from the sun. This 
was accomplished by assuming each flare particle 
independently spirals around the solar magnetic 
field (carried frozen«in by the solar wind), 
while maintaining its pitch angle constant. Among 
the consequences of this transport model are: 

(a) The majority of the solar flare particles 
have pitch angles close to 90°, Calcu- 
lated pitch angle particle distributions 
as a function of time agree fairly well 
with observations at 1 AU. 

(b) The solar flare particle flux rate is 
ndent on distance from the sun 

while the time integrated 
particle flux obeys the expected inverse 
square distance relationship (R"^). 

(c) The solar flare particle energy spectra 
are essentially independent of distance 
from the sun, since particle losses after 
acceleration are neglected. Some parti- 
cle losses probably occur, leading to a 
softening of the energy spectra and 
larger flux spatial dependences than 
herein calculated. 


quite depe 


q Particle charge (coulomb) 

R Distance from sun (AU) 

r Radius of particle gyro rotation around 

solar magnetic field line (meters) 

s Length of solar magnetic field flux line 

(meters) 

t Time (hrs) 

^rise Onset + rise time (time from flare on sun 
to particle fliox maximum observed at R AU) 

( seconds ) 

tj. Particle transit time from the sun (seconds) 

V Particle velocity (meters/sec) 

v^l Component of particle velocity parallel to 

solar magnetic field (meters /sec) 

v^ Component of particle velocity perpendi- 
cular to solar magnetic field (meters /sec) 

O' Parameter in solar flare particle event 
model 


This solar flare particle event model is directly 
applicable to missions which involve sending a 
spacecraft on missions in the ecliptic plane far 
from 1 AU from the sun. It is particularly appli- 
cable to missions to Mercury and to the outer 
planets. It is not intended to apply far from 
the ecliptic plane or past the heliosphere. With- 
in these limits it provides an easy-to-use improve- 
ment over spatially independent models. 

VI . I^ome nc la ture 

A Normalization factor for solar flare 

particle event model 

AU Astronomical unit (l.5 x 10^ km) 

B Solar magnetic field strength (webers/ 

meter^) 

Bj. Radial component of solar magnetic field 
strength (webers /meter^ ) 


H Pitch angle of particle velocity with 

respect to solar magnetic field (degrees) 

4> Particle (proton) flux (particles/cm^-hr- 

Mev) 

/4^ Time-integrated particle flux (particles/ 

cm^) 

A 

^ Peak particle flux rate (part ides /cm^-hr) 

6 Azimuthal angle (polar coordinate) 
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Table T . Annual Totals for golar Flare Radiation Event Particle Fluxes at 1 AU 



Table i. Hniiuaj. xvjua-i-o j. 

Niimbar of Events (J<j>>10 Mev) 

Annual Totals ( particles /cm^) 

Year 

>10^ >10'^ >10® >10^ >10^° 

J4>>10 Mev f<l>>30 Mev /4>>100 Mev 

1956 

11 1 

2.T X 10^ 1.1 X 10^ 3»T X 10® 

1957 

1 4 6 3 

1.0 X 10^° 8.5 X 10® 5.9 X lo"^ 

1958 

1 4 2 

1.1 X 10^° 1.4 X 10^ 9.6 X lo"^ 

1959 

11 5 

2.2 X 10^° 4.3 X 109 4.2 X 10® 

i960 

3 11 1 1 1 

1.3 X 10^° 2.7 X 10^ 3.9 X 10® 

1961 

15 3 1 

2.1 X 10^ ^.2 X 10® k,9 X 10"^ 

1962 


1.6 X 10^ 2 X 10^ 1.4 X 10^ 

1963 

111 

5.6 X 10® 6.5 X lo"^ 6.4 X 10^ 

1964 

.1 1 

1.2 X 10^ 1,2 X lo"^ 1.1 X 10^ 

i 1965 

2 1 

6 

2.2 X 10® 2.2 X 10' 2.0 x 10 

1966 

3 2 

8.1 X 10® 1.0 X 10® 1.0 X lo"^ 

• 1967 

3^11 

2.2 X 10^ 2.T X 10® 2.8 X lo"^ 

1968 

1 3 1 

2.1 X 10^ 2.4 X 10® 2.3 X lo"^ 


Table II. Transit Times (Seconds) for Protons as a Function of Pitch Angle at 1 AU 


i^roton 

Pitch Angle 



Energy 

(Mev) 

p= 0° 

P= 15° 

P = 30° 

P= 45° 

P= 60° 

P= 75° 

P= 80° 

1.0 

1.23 X 10**- 

1.27 X 10^ 

1.42 X 10^ 

1.T4 X 10^^ 

2.46 X 10^ 

k 

4.75 X 10 

7.08 X 10*^ 

3.0 

T.l X 10® 

7.32 X 10® 

8.2 X 10® 

1.0 X 10“+ 

1.42 X 10^ 

2.74 X lO**’ 

4.07 X lo‘*’ 

10 

3.9 X 103 

4.04 X 10^ 

4.51 X 103 

5.52 X 103 

7.82 X 10® 

1.51 X 10*^ 

2.25 X 10*^ 

30 ! 

2.24 X 10® 

2.3 X 10® ! 

2.58 X lo3 

3.16 X 103 

4.47 X 103 

8.62 X 10® 

1.28 X 10*^ 

loo 

1.23 X 10® 

1.27 X 10® 

1.42 X lo3 

1.T4 X 103 

2.46 X 10® 

4.75 X 10® 

7.03 X 10® 

300 

8.18 X 10^ 

8.47 X 10® 

9.45 X 102 

1.16 X 103 

1.64 X lo3 

3.16 X 10® 

4.71 X 10 ® 
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Table III. Minimum Transit Times for Protons as a Function of Energy (p= O) 


Distance 
From Sun 
(AU) 



Proton Enere 

y (Mev) 

1 

2 

10 

30 

100 

1 ^00 

0.1 

1.07 X 10^ 

6.2 

X 

io2 

3.4 X 10^ 

1.95 X 10^ 

1.07 X 10^ 

7.1 3 

clO^ 

0.4 

4.3 X 10^ 

2.48 

X 

10^ 

1.36 X 10^ 

7.8 X 10^ '' 

4.28 X 10^ 

2.84 

X 10^ 

0.7 

8.0 X 103 

4.63 

X 

10^ 

2.55 X 103 

1.46 X lo3 

8.0 X lo3 

5.33 

X 10^ 

i;o 

1.23 X 10^ 

7.1 

X 

103 

3.9 X 103 

2.24 X lo3 

1.23 X 10^ 

8.18 

X 10^ 

( 1-5 

1 

2.09 X 10^ 

1.21 

X 

10^ 

6.6 X lo3 

3.8 X lo3 

2.09 X lo3 

1.39 

X 10^ 

1 2.0 

3.16 X 10^ 

1.82 

X 

10^ 

1.0 X 10^ 

5.74 X lo3 

3.16 X 103 

' 2.1 X 

103 

3.0 

6.05 X 10^ 

3.48 

X 

10^ 

1.92 X 10*^ 

1.1 X 10^ 

6.05 X lo3 

4.0 X 

10^ 

4.0 

9.95 X 10*^ 

5.75 

X 

10^ 

3.16 X 10*^ 

1.81 X 10*'' 

9.9 X lo3 

6.6 X 

103 

5.0 

1.5 X 10^ 

8.65 

X 

10** 

4,75 X lO**- 

2.73 X lO*' 

1.5 X lO*' 

9.95 : 

X 103 


Table IV. Transit Times (Seconds) for 1 Mev Proton as a Function of Pitch Angle 


Distance 
Prom Sun 

/ att\ 



Pitch Angle 

(AU) 

p= 0° 

P= 15° 

p= 30° 

P= 45° 

P= 60° 

P= 75 *^ 

p= 80° 

0.1 

1.07 X 10^ 

1.11 X 10 ^ 

1.24 X 10 ^ 

1.52 X 10 ^ 

2.14 X 10 ^ 

4.13 X 10 ^ 

6.15 X 10^ 

0.4 

4.3 X 10 ^ 

4.45 X 10 ^ 

4.97 X 10 ^ 

6.1 X lo 3 

8.6 X iq 3 

1.66 X 10 ^ 

2.48 X lO*' 

0.7 

8.0 X iq 3 

8.27 X 103 

9.25 X 103 

1.13 X 10*^ 

1.6 X 10 *'' 

3.09 X lO*' 

4.6 X 10 **' 

1.0 

1.23 X 10 ** 

1.27 X 10 ^ 

1.42 X lO*' 
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Figure 1. Characteristics of the Quiet-Sun 
Magnetic Field Carried Frozen in hy the Solar 
Wind, 



Distance From the Sun (AU) 

Figure 2. Proton Gyro- Radius as a Function of 
Distance From the Sun for Undisturbed Inter- 
planetary Magnetic Field. 



and Proton Energy 








AU Values 










